When gel and glass meet: A mechanism for multistep relaxation 
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We use computer simulations to study the dynamics of a physical gel at high densities where 
gelation and the glass transition interfere. We report and provide detailed physical understanding 
of complex relaxation patterns for time correlation functions which generically decay in a three-step 
process. For certain combinations of parameters we find logarithmic decays of the correlators and 
subdiffusive particle motion. 
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Various types of phase transformations in disordered 
materials have been described in the literature, such as 
percolation, gelation, or glass transition and it is 
thus natural to seek novel behaviors that emerge from 
the competition between these different transformations. 
Recently, the structure and dynamics of mixtures of large 
and small colloidal particles [2|, y] , of materials with com- 
peting lengthscales 0, Q , or of systems with competing 
interactions 0, 0] has been studied using a variety of the- 
oretical and experimental techniques. The motivation to 
carry out these investigations is of fundamental as well 
as practical nature, since new phenomena can been dis- 
covered and disordered materials with novel properties 
are created. 

Of particular interest are the similarities, differences 
and competition between gelation and glass transition 
which can be observed in a number of soft materi- 
als [1, 0, [n], Previous work focused on particle 
systems with a hard core repulsion competing with a 
very short-range attraction [6( , which can be realized ex- 
perimentally in colloid-polymer mixtures Such sys- 
tems can be handled theoretically using liquid state the- 
ory for the structure, and mode-coupling theory for the 
dynamics, which then yields detailed predictions for the 
location and nature of the dynamic transitions @, [l2j ]. 
Non-ergodic solids are predicted both upon compres- 
sion (glass physics) or upon increasing the attraction 
strength (gelation), with a peculiar dynamical behav- 
ior (logarithmic decay of time correlation functions and 
sub-diffusive particle displacements) in the region of pa- 
rameter space where these transition lines intersect. A 
number of numerical simulations 
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., UM an d experi- 
ments d, US U) Hi have confirmed the overall topology 
of the phase diagram, with some simulations reporting 
the dynamical signatures predicted theoretically jH, [3] ■ 
However, a detailed physical understanding of these dy- 
namics and an exploration of their universality for more 
complex materials has not been achieved. Moreover, 
since mode-coupling singularities are generically avoided 
in real materials |l|, it is important to explore also the 



robustness or generality of these results [16| , and the pos- 
sible deviations or new processes which might emerge in 
real materials. 

We explore these important open directions using com- 
puter simulations of a model system that is a coarse- 
grained representation [l7| of a transient gel which has 
been studied in experiments [l8| . In this system an equi- 
librium low-density gel is obtained by adding telechelic 
polymers to an oil-in-water microemulsion. Since the 
polymer end-groups are hydrophobic, the polymers effec- 
tively act as (attractive) bridges between the oil droplets 
they connect, whose strength, lengthscale and typical 
lifetime can be controlled at will. Denoting by the 
number of polymers connecting droplets i and j we have 



established in Refs. [Tfl, LL9| that the following interac- 



tion is an effective coarse-grained representation of this 
ternary system: 



V = J2 (—) +ei5Z<^V FENE M+eo^a i . (!) 
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The first term is a soft repulsion acting between bare oil 
droplets, where cry = (a - ; + <jj)/2, <Ji is the diameter of 
droplet i, and ry is the distance between the droplet cen- 
ters. The second term describes for the entropic attrac- 
tion induced by the telechelic polymers, which has the 
standard "FENE" form known from polymer physics [2(| , 
Vfene(^) = hi(l — (r,y — <7ij) 2 /£ 2 ), and accounts for the 
maximal extension I of the polymers. The last term in- 
troduces the energy penalty eo for polymers that have 
both end- groups in the same droplet. The most drastic 
approximation of the model (P) is the description of the 
polymers as effective bonds between the droplets, which 
is justified whenever the typical lifetime of the bonds is 
much larger than the timescale for polymer dynamics in 
the solvent 18j . In order to describe the dynamics of the 
system, we use a combination of molecular dynamics to 
propagate the droplets with the interaction |T]), and lo- 
cal Monte Carlo moves with Metropolis acceptance rates 
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FIG. 1: Time dependence of F s (q,t) for q — 6. a) Gelation 
at <j> = 0.5, Tii n k = 10 2 and increasing R. b) Glass transition 
at R = and increasing cj>. 

Tjj^ k min[l, exp(AV/kBT)] to update the polymer connec- 
tivity matrix (ItI fl9j | . Thus Tn n k is the timescale 
governing the renewal of the polymer network topology. 
In order to prevent crystallization we use a polydisperse 
emulsion with a flat distribution of particle sizes in the 
range en € [0.75, 1.25]. 

For moderate volume fraction, (j> ss 0.2, the model be- 
haves as an equilibrium transient gel with non-trivial dy- 
namics [13] if R, the fraction of polymers per oil droplet, 
is larger than the percolation threshold R p « 2. In this 
case, the system is viscoelastic with a relaxation time set 
by the bond lifetime Ti; n k, and mechanical strength set 
by R, as illustrated in Fig. [T^ where we show the time 
dependence of the self-intermediate scattering function, 

F * («> *) = jf E ex P ^ ■ (*) ~ v i (°)1 ) > ( 2 ) 
j 

across percolation for <^ = 0.5 (q = 6.0, near the main 
peak in the static structure factor). A two-step decay is 
observed with a plateau height controlled by R, and a 
slow decay at times t s» Tu n k. These features reflect the 
vibrations of an increasingly stiffer network of connected 
particles, followed by a slow reorganization of the tran- 
sient network. In the opposite limit where R — and the 
volume fraction (j> becomes large, the microemulsion be- 
comes a standard dense glass, Fig.QJi. Here, the two-step 
decay stems from particle vibrations within the transient 
cage formed by the neighbors, followed by slow structural 
relaxation. As usual in this situation, the relaxation time 
increases dramatically with <p while the plateau height re- 
mains constant 

Our goal in this work is to explore the competition 
between the two well-documented phenomena illustrated 
in Fig. fn The space of control parameters is large, so 
we fix 17] {£ = 3.5, k B T = e = l,ei = 50} and vary 
{(f), R, Tr m k}. We successively describe the effect of (i) 
increasing the density of an equilibrium gel, (ii) adding 
attractive interactions to a viscous close to the glass tran- 
sition, (iii) changing the bond lifetime of a system close 
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FIG. 2: Compressing an equilibrium transient gel with R — 8 
towards the glass transition. The symbols are for rii n k = 10 2 
and the full lines are for frozen bonds, nink = oo. The dy- 
namical evolution is a non-trivial mixture of both gelation and 
glass transition with a logarithmic decay towards a plateau at 
intermediate volume fraction. 

to dynamic arrest to estimate the relative importance of 
bonding and steric hindrance for the relaxation. 

We start our exploration in Fig. [2] which shows the 
evolution of the equilibrium gel dynamics, R — 8, rii n k = 
10 2 , as it is compressed from <p — 0.5 towards the glass 
transition. Contrary to the extreme cases shown in Fig. [2 
here both the relaxation time and the plateau height 
increase simultaneously. Thus these results are not ex- 
plained by gelation or the glass transition alone, but truly 
result from the non-trivial effect of their competition. 
This is directly demonstrated by increasing, for the same 
state points, the bond lifetime Tr ln k to a very large value. 
For = 0.5, the correlator quickly decays to the plateau 
and then becomes completely arrested, showing that at 
low density only the transient network physics is at play. 
At larger density, <j> — 0.60 and 0.65, a very slow decay 
towards a plateau is observed. In particular, for 4> = 0.65 
the data follow a nearly logarithmic decay over about 5 
decades in time, a behavior which is seen neither for the 
gel nor the glass alone. This means that the elastic re- 
laxation of the network is considerably slowed down by 
crowding effects. Interestingly, qualitatively similar ex- 
perimental observations have recently been reported in 
micelles [llj and attractive nanoparticles [l5|. At even 
larger volume fractions, <f> = 0.7, F s (q,t) decays again in 
a two-step process, but the slow decay is now controlled 
both by the bond-lifetime (as in gels) and by density (as 
glassy liquids). Although the dynamics at large 4> and 
R is a two-step process similar to the one of the glassy 
fluid at R = 0, the effect of Ti in k on the second decay es- 
tablishes its very different nature. While crowding alone 
is responsible for the slow dynamics near R = 0, bond- 
lifetime and network reorganization control the dynamics 
at large R. This physical distinction is reminiscent of the 
"bonded glass" and "repulsive glass" nomenclature intro- 
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FIG. 3: Effect of increasing the attractive interactions, R, in 
a dense repulsive fluid, <f> = 0.61, for a finite bond lifetime, 
Tiink = 10 2 on (a) F s (q,t) and (b) A 2 (t). For intermediate R 
values, a nearly logarithmic time decay of F a (q,t) is observed 
(symbols), associated with a sub-diffusive behavior of single 
particle displacements. 



duced for attractive colloids [2l| . It is between these two 
regimes that multistep relaxations can be observed. 

Alternatively, gel and glass can meet by increasing the 
attraction in dens e sy stems, as usually done in polymer- 
colloid mixtures [lfj, [l3|. Note that in our model we 
can increase the attraction between droplets, R, without 
simultaneously increasing the bond lifetime Tn n k, while 
these quantities are tightly coupled in attractive col- 
loids [TO, EH, HH]. This allows us to disentangle static 
from dynamic effects. In Fig. [3] we show the effect of in- 
creasing R in a dense fluid at (j> = 0.61 using a constant 
bond lifetime rank = 10 2 . We show both the evolution of 
F s (q,t), and of the mean-squared displacement, 



A 2 (i) 



^X>«(t)-r<(0)| a ) 



(3) 



As for the low density system in Fig. [TJ increasing the 
attractive interaction allows the system to cross the per- 
colation line and to become viscoelastic, i.e. F s (q,t) de- 
cays in two steps. However, the proximity of the glass 
transition makes the dynamic evolution more complex, as 
this produces, at intermediate values of R « 5 — 7 a very 
slow decay of the time correlation function, which can 
be empirically described by a logarithmic decay, Fig. [3^. 
This behavior resembles the ones reported in numerical 
work for colloids with short-range attraction 13j, Il4j . 



The mean-squared displacement also evidences devia- 
tions from a simple two-step process, as particle displace- 
ments appear to be transiently sub-diffusive at those in- 
termediate R where logarithmic decay is observed. Again 
this behavior can not be accounted for by a simple "su- 
perposition" of gel and glass dynamics. 

We explore further the dynamics at those state points 
where the dynamics is non-trivial, and make use of the 
flexibility offered by the present model to change attrac- 
tion and bond-lifetime independently. In Fig. [H we show 
the effect of increasing ru n k on F s (q, t) and A 2 (i) for con- 




FIG. 4: Increasing the bond lifetime runk from 1 to 10 4 at <f> = 
0.61 and R = 6.3 produces a crossover from nearly logarithmic 
(dashed line) to a three-step process when timescales are well- 
separated. The inset shows the mean-squared displacement 
for the same parameters, revealing the existence of the distinct 
lengthscales (dashed lines) controlling the dynamics. 

stant values of R — 6.3 and tfi — 0.61, i.e. where the log- 
arithmic decay was most prominent in Fig. [3j We now 
realize that the logarithmic decay seen for rii n k = 10 2 
is in fact a very specific instance of a more generic 
three-step decay of time correlation functions which re- 
flects the double localization of the particles within their 
cages and within the transient particle network. This 
three-step process is also clear from the behavior of the 
mean-squared displacements (inset of Fig. 2]), which re- 
veals the existence of two lengthscales controlling the dy- 
namics of this system. A first plateau is observed for 
A 2 k 0.1, which corresponds to the typical cage size in re- 
pulsive glasses [l[ , and a second plateau is observed near 
A 2 w 0.5, which corresponds to the localization of the 
particle due to the presence of the percolating polymer 
network 17j. The logarithmic decay is only observed if 
a non-generic combination of timescales and lengthscales 
combines both slow processes in a single, nearly logarith- 
mic one. In fact our data suggest that the generic situa- 
tion should be the occurrence of a three-step process cor- 
responding to well-separated bond/cage relaxation. To 
the best of our knowledge, there is no experimental report 
of such a three-step decay of time correlation function. 
We s ugg est that compressing the transient gel system of 
Ref. [l8( is a possible route to such observations. 

We now summarize our findings in a phase diagram in 
Fig. El The low density part, < 0.2 was described ear- 
lier [171 ]. It contains a sol phase if <j> and R are small, and 
a phase separated region when attraction is increased at 
low <j). There is a large region (currently explored in ex- 
periments 0, H3|) where a homogeneous, equilibrium, 
transient gel is formed. In the present work we have ex- 
plored the interplay between the glass transition at large 
volume fraction and gelation at large R, where we dis- 
covered new dynamical phenomena. The iso-relaxation 
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gations. Even more complex relaxation dynamics could 
potentially emerge in transient gels made with polymer 
mixtures having distinct lengthscales, thus opening the 
door to the creation of new materials. 

Financial support from Spanish project FIS2009-08451 
and ANR's TSANET and SYSCOM is acknowledged. 



FIG. 5: Phase diagram for rji n k = 1.0: Coexistence region 
(o), sol (O) an d gel (□) phases, and measured glass line (A), 
determined from a Vogel-Fulcher fit of the relaxation time. 
Also shown are the iso-r lines (v) for r = 10, 30, 10 2 , • • ■ , 10 4 . 



time lines reported in Fig.[5j obtained for rii n k = 1.0 and 
defined by F s (q, r) = 0.03, show that the dynamics slows 
down both by increasing <j> or R. In the region where both 
R and <j> act as a source for slow dynamics we generically 
obtain a three-step relaxation process, which produces 
for specific combinations of R, <fi, and TH n k, a logarithmic 
decay towards a plateau followed by a slower decay, as 
in Fig. [21 or a fully logarithmic decay as in Fig. [3] The 
present model thus captures a broad range of behaviors, 
and also predicts new types of relaxations. 

The phase diagram in Fig. [5] shows that, in the gel 
phase, ergodic behavior is found even if the attraction 
gets very large. For (f> = 0.5, we have done simulations 
up to R — 50 and still found ergodic behavior, the re- 
laxation time increasing smoothly with R. Thus we find 
no evidence of an 'ideal' gel phase, as predicted theo- 
retically for attractive particles [H, [13, Hj|, although of 
course the relaxation time can get very large if both R 
and Tii n k increase. This suggests a fundamental difference 
between gelation and glass transition because in (fragile) 
glass-forming materials, the predicted mode-coupling sin- 
gularity is also avoided, but it is believed to be replaced 
by a truly divergent timescale [lj . As suggested also from 
simulations of particles with attractive patches 24j, gels 
could be the analog of Arrhenius (strong) glasses with no 
finite temperature singularity. 

One important finding of our study is the prediction 
of a generic three-step decay of time correlation func- 
tions, which has not yet been directly observed experi- 
mentally. As shown in Fig. [4J this requires the existence 
of two distinct lengthscales and well-separated timescales 
for cage/bond relaxation. We suggest that the obser- 
vation of a two-step process for yielding in a sheared 
colloid-polymer mixture [25| might be a rheological ana- 
log of Fig. [H which certainly calls for further investi- 
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